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I. INTRODUCTION 

The single electron transistor is a widely used de- 
vice technology for the exploitation of quantum mechani- 
cal effects (Coulomb blockade) for applications including 
metrology^' and quantum information processing.!^ The 
first devices were fabricated using a combination of a 
metal (aluminium, to form the leads and the island) and 
an insulator (aluminium oxide, to form the tunnel bar- 
riers) P*^ such devices are known as metallic single elec- 
tron transistors (metal SETs). Much of the understand- 
ing of single electron effects in this type of device was 
obtained by measuring their electrical characteristics as 
a function of temperature and magnetic field. Today, 
metal SETs have a well understood behaviour, exhibit- 
ing regular Coulomb oscillations (constant gate period 
and peak height) and a temperature dependence that is 
consistent with the charging energy given by the experi- 
mental Coulomb gap. 

The difficulty in fabricating these lithographically de- 
fined metal island SETs has largely limited their opera- 
tion to very low temperatures (usually requiring a dilu- 
tion refrigerator), so that many other material combina- 
tions have been explored to find an alternative fabrication 
route. One of the most successful has been the degener- 
ately doped silicon single electron device, fabricated us- 
ing silicon-on-insulator technology.— This material system 
was used to demonstrate some of the earliest examples 
of single electron classical logicP and memorjC! in semi- 
conductor devices, as well as providing a route towards 
qubit implementation.'Sl This approach offers significant 
scalability and operating temperature advantages com- 
pared to most metal SET systems. 
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In contrast with metal SETs, the electrical character- 
istics of degenerately doped silicon SETs usually exhibit 
a disordered characteristic (varying gate period and peak 
height). Furthermore, single electron transistor charac- 
teristics are observed without the need for a different ma- 
terial to form the tunnel barriers, as all of the island, the 
leads and the tunnel barriers are made from degenerately 
doped silicon. In semiconductor materials a very effective 
barrier can be formed by free carrier depletion, however 
such a depletion process is significantly constrained by 
the very high doping density used to fabricate these de- 
vices. The doping densities used for degenerately doped 
silicon SET fabrication are very much higher than situa- 
tions where depletion is usually employed. The minimum 
doping level is set by the need to be above the Mott- 
insulator transition (« 5x lO-'^^cm"'^ for P doped silicon) 
to avoid carrier freeze-out at low temperature and the 
maximum doping level is set by the solid solubility limit 
(w lO^^cm"^ for P doped silicon). Coulomb blockade de- 
vices have been successfully dem onstrated over the whole 
of this doping density rang^Sllfll -^^jth lateral dimensions 
differing only by a factor of ~3. It is interesting that 
a change in doping density of more than two orders of 
magnitude makes little difference in the size or geometry 
of devices needed to show Coulomb blockade character- 
istics. 

The earliest silicon SETs were fabricated with con- 
strictions to defi ne the location and strength of the tun- 
nel barriers,'^'^ but it was soon found that these were 
un-necessary for the purpose of demonstrating Coulomb 
blockade characteristics!-^ At that time, it was thought 
that fluctuations in the surface potential on the side walls 
of the silicon were the origin of the tunnel barriers. This 
idea was supported by the change in characteristics ob- 
served each time the device was cycled between the room 
ambient and the measurement temperatures, on the prin- 
ciple that the surface traps would be refilled in some ran- 
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dom way each time. Further support for this idea comes 
from the observation that the temporary apphcation of 
a large bias voltage to the substrate (particularly during 
cool-down) is effective in shifting the overall characteris- 
tics of the system, say to decrease conduction in a device 
where the Coulomb peaks are poorly defined. 

The disorder in the Coulomb peaks makes it difficult 
to extract an unambiguous charging energy. Devices of 
this type can show a very wide range of charging energies, 
within the same transistor or within a batch of similar 
devices, even when fabricated together. Furthermore, in 
contrast to metal SETs, their temperature dependence 
is not well predicted from the charging energy and stan- 
dard theory. Indeed, devices made in this laboratory have 
experimentally demonstrated maximum operating tem- 
peratures from a few tens of Kelvin up to nearly room 
temperature, despite being fabricated to the same speci- 
fications and from the same material. 

Finally, the more recent observation of high quality fac- 
tor resonances during broadband microwave spectroscopy 
is another dramatic difference between the behaviour of 
degenerately doped silicon SETs and metal SETs. Previ- 
ous measurements on metal SETs showed only low qual- 
ity factor resonances (quality factor < 100) due to stand- 
ing waves within the waveguide system used to couple 
the microwave source to the device under test.^^^ These 
standing waves cause the amplitude of the microwave 
signal at the device to vary with frequency and the 
device current is influenced by this amplitude through 
heating effects. The current passed by a device under 
fixed source-drain and gate bias increases with temper- 
ature in a characteristic way due to the co-tunnelling 
effect. This is confirmed by the experimental obser- 
vation in the case of metal SETs that the device cur- 
rent due to the microwaves always increases and that 
the frequency dependence is unaffected by the bias con- 
ditions of the device under test. This is in contrast 
to the experimental observations on the degenerately 
doped silicon SETs, where both increases and decreases 
in device current due to the microwaves were seen and 
the frequency dependence was influenced by the device 
bias conditions.li^ The quality factors of these resonances 
can be very high (up to around 400,000) and the peak 
shape varies from positive/negative gaussian to differen- 
tial. The very large number of high quality factor res- 
onances (more than 1000 seen in the frequency range 
l-20GHz) and the Rabi-like oscillations under pulsed ex- 
citation conditions.'^ strongly indicates a quantum me- 
chanical origin for these effects. 

In this paper, we attempt to explain the behaviour of 
degenerately doped silicon SETs on the basis of local- 
ized transport resulting from significant charge trapping 
on the sidewalls during cool-down. This charge trapping 
is influenced by the details of the device fabrication and 
can explain the relative independence of the required de- 
vice geometry on the doping density. We will firstly de- 
scribe the observed temperature dependence and disor- 
dered Coulomb blockade response at DC in terms of hop- 
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FIG. 1: Schematic diagram of the radiation coupling method 
using separate DC and AC lines. The DC lines are low pass 
filtered (LPF) by MiniCircuits BLP1.9 units. Top inset: SEM 
image of a typical SET. Bottom inset: schematic cross section 
of the silicon-on-insulator structure. 



ping transport via impurity states. Next, we will discuss 
how trapped charges can give rise to high quality fac- 
tor resonances during microwave spectroscopy and what 
characteristics these resonances might be expected to ex- 
hibit. 



II. EXPERIMENTAL DETAILS 

Devices investigated in this work were fabricated us- 
ing degenerately P-doped silicon-on-insulator wafers (see 
bottom inset of Fig. [T]). Devices were fabricated using 
electron beam lithography and a sacrificial metal mask 
to transfer the pattern by reactive ion etching. Patterned 
substrates were thermally oxidized at 1000°C to reduce 
the surface trap state density, which in turn reduces the 
incidence of random telegraph signals. A micrograph of 
a typical measured SET is depicted in the top inset of 
Fig.[l] 

Measurements were carried out in a specially con- 
structed probe that could be directly immersed in a liquid 
cryogen, so that the device under test and a substantial 
part of the connecting leads, were all efficiently cooled. 
All measurements discussed here were performed in liquid 
helium, at a temperature of 4.2K. DC measurements were 
made using Keithley 236 source-measure units (SMU) 
controlled by a Lab View programme via the GPIB bus. 
The DC measurement connecting leads were directly cou- 
pled to the source, drain and gate terminals of the device 
and individually low pass filtered by MiniCircuits BLP1.9 
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inline units at room temperature with cut-off frequency 
of 1.8 MHz. This filtering is essential to prevent elec- 
tromagnetic noise (due to external sources) from propa- 
gating to the device under test, which would otherwise 
be subject to strong electron heating. Microwave signals 
were provided by an Agilent E8257D source capable of 
delivering up to -)-15dBm via a flexible sma-terminated 
cable to the probe. 

AC signals can be applied to the source-drain leads 
via a biasing circuit; a schematic of the coupling circuit is 
shown in Fig.[T] A 50 metal film resistor is provided to 
impedance match the waveguide and microwave source, 
so that standing waves are inhibited. 100 uF ceramic 
capacitors are used to decouple any DC voltage on the 
waveguide from reaching the device. At 100 MHz (the 
lowest frequency used here) the impedance of each ca- 
pacitor is < 1 and decreases with increasing frequency 
so that the voltage drop across the 50^1 load resistor ap- 
pears as a voltage difference between the source and the 
drain sample connecting leads, with negligible loss. The 
10 /iH inductor (impedance > 6 kfi at 100 MHz) prevents 
the AC voltage from reaching the SMU, where it could 
interfere with the DC current measurement. 

Investigations have been previously carried out where 
the radiation was indirectly coupled to the sainple by 
broadcasting it from an open-ended wave-guide.'^ This 
indirect coupling method resulted in a global time depen- 
dent potential for the whole device, whereas the more ef- 
ficient coupling technique presented here results in a time 
dependent potential difference between device electrodes 
(mainly source and drain). 



III. DC TRANSPORT 

A typical stability plot, measured at a temperature of 
4K, is shown in the inset of Fig. [2jb) showing aperiodic 
Coulomb oscillations and unequal peak heights. Clearly 
this device must meet the two requirements for the ob- 
servation of Coulomb blockade; the thermal energy must 
be smaller than the charging energy and the tunnel bar- 
rier impedance must be large compared to the resistance 
quantum. The second condition requires that the tun- 
nel barrier impedance should be at least lOOkfi and the 
first condition requires that the tunnel barrier does not 
extend by more than one or two squares along the de- 
vice axis (otherwise the capacitance of the island is too 
large). However, a doping density of falO^® cm^'^ gives 
a sheet resistance of about IGOOrJ/D at room tempera- 
ture, dropping by a factor of about 1.6 at 4K^ giving 
a "tunnel barrier" resistance of only a few kil from geo- 
metrical considerations alone. Clearly, the carrier trans- 
port in the tunnel barrier region must be influenced by 
some other mechanism that greatly increased the resis- 
tance. To explain the magnitude of the current observed 
at a Coulomb peak, the tunnel barrier resistance must 
be nearly two orders of magnitude larger than the value 
expected from the geometry and from the doping density. 



It was previously suggested that the high tunnel barrier 
resistance is a result of the depletion of free carriers by 
charge trapped at the sidewalls of the silicon.!^ Of course, 
a high density of traps at the sidewalls is expected due to 
the non-(lOO) orientation of the surface.^* In general, the 
density of traps due to unsatisfied dangling bonds at the 
silicon surface is dependent on the crystalline orientation 
of the sidewalls. Oxidation is routinely used to passi- 
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FIG. 2: Depletion effect due to sidewalls trapping, a) Cartoon 
showing a top-view of the nano-structure. Depleted regions 
are mostly formed in the vicinity of the side edges and where 
the channel is narrower, b) Cross sectional electrostatic po- 
tential at the generic position marked by the dotted line in 
(a). Fermi energy is bent because of charge accumulation at 
the sidewalls. Ionised donors and neutral sites with bound 
electrons are represented by plusses and crosses, respectively. 
Inset: a device stability plot: drain-source current as a func- 
tion of drain-source voltage, Yds, and side gate voltage, Ya- 
Blue-to-red scale corresponds to current varying from -60 pA 
to -1-60 pA. 
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vate such dangling bonds, but little difference in device 
behavior was seen for the earliest devices made without 
oxidation and those with oxidation. Oxidation is effec- 
tive in reducing random telegraph noise, but has little 
effect on device geometry necessary for the observation 
of Coulomb Blockade. Moreover, the depletion depth due 
to these traps would be fixed (and estimated to be much 
too small to explain the tunnel barrier size) for a given 
doping density. So that an increase in the doping density 
would require a corresponding decrease in the physical 
size of the tunnel barrier, contrary to experimental ob- 
servation. Another mechanism for trap formation arises 
from the fate of dopant atoms in the silicon etched away 
to form the sidewalls. Reactive ion etching is used to form 
the island, tunnel barriers and leads by removing the un- 
wanted silicon. This process proceeds by transporting 
gas reactants to the silicon surface, where a chemical re- 
action takes place to form volatile products, which are 
removed by the vacuum pumping system. Typical sil- 
icon dopants, such as phosphorus, do not form volatile 
compounds with the CF4 and SiCU reactive etch gases 
at the usual etch temperatures. The phosphorus atoms 
remain on the surface after etching and become incor- 
porated into the oxide during passivation. These atoms 
are no longer able to act as dopants, but do act as traps. 
The density of these traps directly depends on the dop- 
ing density of the silicon so that a change in this doping 
density has no direct effect on the depletion depth in the 
silicon. This mechanism explains why all tunnel barriers 
appear similar for a wide range of doping density. 

Transport in highly doped silicon SETs is strongly lo- 
calized in and around the SET island and tunnel barri- 
ers, due to depletion resulting from trapped charge on 
the sidewalls. At high temperatures, these traps will be 
mostly empty due to thermal excitation; so that the re- 
sistance is little changed from that expected on the basis 
of the bulk material properties. But, as the temperature 
is reduced these empty, neutral traps are able to capture 
electrons that tunnel in from the nearby silicon, which 
results in a significant negative charge accumulation at 
the sidewalk To maintain overall charge neutrality, this 
negative charge causes depletion of free carriers from the 
nearby silicon; these being the most likely to be captured 
in the traps anyway. This depletion increases the size 
of the tunnel barrier; but in practice, would only affect 
the resistance of narrow sections, such as the SET island 
and tunnel barriers (Fig. [2]) . Such depletion moves the 
Fermi energy down into the impurity band, despite the 
high doping density, so that the conduction mechanism 
changes from metallic to a hopping behaviour. 

Such conduction behaviour normally arises only in 
lightly doped semiconductors at low temperatures when 
the previously free carriers are recaptured by the ionised 
dopants and the overlap between carrier wavefunctions 
on adjacent impurity sites becomes weak. But such char- 
acteristics can also be observed in systems where the car- 
rier density is usually high but is subject to depletion 
due to a gate voltage or by trapped charge. In that case. 



the conductivity decreases exponentially as the energy or 
spatial separation between the localized sites increases. 
We can apply the analysis of Efros and ShklovskiP^l to 
these situations. 

Fig. |3] shows the electrical differential conductance 
measured as a function of temperature and source-drain 
voltage for a highly doped silicon nano-wire device, un- 
der conditions where the Coulomb Blockade effects will 
be mostly suppressed by the applied potential difference 
at drain-source and/or gate. Measurements have been 
performed at six different temperatures for each applied 
electric field by means of standard lock-in techniques. 
The data points shown are the average values for the 
differential conductance from a number of equivalent ex- 
perimental conditions obtained by sweeping Vds in both 
forward and reverse directions. The observed behaviour 
is typical of a wide range of device types, including those 
with lithographically defined tunnel barriers. It is clear 
that both temperature and electric field have a strong 
influence on the transport properties. This current de- 
pendence is often attributed to co-tunnelling in Coulomb 
Blockade devices, where the tunnel barrier is formed by 
a material with a large energy gap. An analytic form 
for the inelastic co-tunnelling contribution to the sample 
conductivity of an SET measured at low temperature T 
and applied source-drain voltage Yds has been given by 
Averin and NazaroA^ as 

acot - A{BT^ + CVEs) (1) 

where A, B and C are constants that depend on the na- 
ture of the tunnel barriers and fundamental constants. 
However, it is clear that such a relation does not satisfac- 
torily explain the dependencies shown, where the largest 
changes occur at the lower voltages and temperatures. 

As discussed earlier, the electron density close to the 
sidewalls is expected to be significantly reduced by deple- 
tion resulting from electrons trapped in the grown oxide, 
so that it is appropriate to consider the temperature and 
electric field dependence of the hopping conductivity in 
low-doped semiconductors to explain the characteristics 
in Fig. |3] The conductivity of a disordered system in 
the smaU polaron hopping regime has been discussed by 
Triberis.'^IlThis work concludes that Mott's law gives the 
main temperature dependence and that the electric field 
dependence is given by a function of the square of the 
source-drain voltage. Using these relations, the following 
function has been fitted to the data in Fig. [3] 

ln[a{T,VDs)]^A + BT-^/^ (2) 

where A and B are treated as different fitting constants 
for each value of Vds- The functions that model the field 
dependence are then fitted as follows 

A = C + D{E + VDsr^ (3) 
B = FiE + Vosy^ (4) 
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being C = -1.06±0.08, D = 0.024±0.001, F = 
-0.067±0.001 and E = 0.054 fitting parameters {E is 
used as a fixed parameter, therefore it does not come with 
an error). The fit for low temperatures and small source- 
drain voltages is very satisfactory; at higher voltages and 
temperatures the hopping model becomes progressively 
less applicable. This agreement between the experimen- 
tal data and the hopping model suggests that the ob- 
served temperature dependence is mostly determined by 
changes in the transparency of the tunnel barriers rather 
than by co-tunnelling. 

Very large charging energies have been estimated in the 
past from electrical measurements at low temperatures 
suggesting that operation at high temperature should be 
possible.^' However, these estimates often disagree with 
the observed maximum operating temperature. In these 
cases, it is the degradation of the tunnel barrier resis- 
tance that is the cause of the reduction in maximum op- 
erating temperature, rather than the onset of significant 
co-tunnelling. Longer constrictions seem to give higher 
maximum operating temperatures than shorter or nar- 
rower ones. This effect may be explained by the increase 
in the number of hops taken to pass through the tunnel 
barrier. An increase in the number of hops will reduce 
the energy barrier per hop and so reduce the effect of an 
increase in temperature or electric field. 

In a metallic SET, the electrical characteristics are de- 
termined by the island capacitance and the tunnel barrier 
transmission. The charging energy for an extra electron 
to occupy the island has only one contribution, that due 
to the island capacitance. In a degenerately doped silicon 
SET, hopping transport can modify this behaviour, since 
the energy for an electron to occupy a localized site in 
the island now has three contributions. The first is the 
conventional Coulomb charging energy due to an amount 
of charge Q on an island capacitance C (we will assume 
a fixed capacitance), where the localized site is situated. 
The second contribution is due to the binding energy of 
the site, which will depend on the nature of the localisa- 
tion. This contribution decreases the system energy by 
an amount for each electron. The binding energy for 
an isolated phosphorus atom in a silicon lattice is given 
by the Bohr model as « 44meV. However, in degenerately 
doped material the donor spacing is sufficiently close that 
wavefunction overlap strongly influences this result. Dis- 
order in the donor lattice position leads to a deviation 
from the average donor separation, r, that varies with 
r^] so that while most donors are separated by a spac- 
ing that is close to the average value, a small fraction 
have significantly smaller spacings. At these locations 
the binding energy for one electron is increased and for 
the second is decreased. The effects of the disorder in 
the donor location and the local influence of the trapped 
charge causes the binding energy to be different at each 
localized site (Fig. [2]) . The third contribution is due to 
the electron interactions between the occupied localized 
sites. If the spacing between an occupied site pair is r, 
then this increases the energy by an amount Ei which 
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FIG. 3: Natural log of a device differential conductance at 
different temperature and electric field. Measurements were 
taken at applied source-drain voltages of 0, 0.02, 0.04, 0.06, 
0.08, 0.1, 0.12 and 0.14V from bottom to top. Solid lines are 
the results from a fit using eq. Inset: total energy on the 
island as a function of the stored charge for a classical metal 
Coulomb Blockade device (dotted line) and a highly doped sil- 
icon device (solid line) . The red dots indicate positions corre- 
sponding to integer charge. Of note is that the solid parabola 
is segmented and displaced due to the additional contribu- 
tions from the binding and interaction energies. Some of the 
integer charge positions have multiple solutions corresponding 
to different charge configurations. 

depends on the electron charge e, on the permittivity of 
the channel material e and on r in the following way. 



The total energy of the system Ejot, then depends on 
these three contributions 

2\{N-2)\ N 2 

Etot= J2 ^«w-E^''^w + ^ (6) 

i=l 1=1 

where N is the number of electrons localized at donor 
sites within the island. For the purposes of assessing the 
Coulomb blockade charging energy, we are only interested 
in the change in the total energy Etot for the addition or 
subtraction of one electron. 

e2 

AEtot^NE,{r)-Eu + — (7) 

So that the classical Coulomb Blockade energy is modi- 
fied by two additional terms. The magnitude of the first 
of these additional terms (the interaction energy) will be 
N times about 3 meV (assuming a separation of about 
half the typical dot diameter of «40 nm) and the second 
(binding energy) term will be of the order of 44 meV from 
the previous discussion. Clearly, with about 15 electrons 
on the island these two additional terms are expected to 
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cancel and the classical Coulomb Blockade charging en- 
ergy depends only on the island capacitance. However, 
as considerable variation from these average approxima- 
tions are expected, then the actual charging energy will 
contain a correction term to the classical result, which is 
due to the difference between the binding energy and in- 
teraction terms. This difference will vary as the electron 
number varies and for a given electron number will also 
vary with configuration i.e. if a different set of donors is 
occupied with the same number of electrons then the en- 
ergy difference is also changed. In addition, the estimate 
for the interaction term makes no allowance for screen- 
ing, which will greatly reduce the influence of charges 
at the larger separations. Indeed, the screened Coulomb 
interaction energy will decay rapidly with distance; we 
describe it as follows: 

Escr — EiC ^'"^ (8) 

where l/kg is the Thomas- Fermi screening length. We 
estimate 1/fcs ~114nm for our system; this may increase 
the number of electrons on the island to achieve a balance 
between the binding energy and interaction energy terms 
to the order of 100. The effect of these additional terms is 
illustrated schematically in the inset of Fig. [3j where the 
parabolic dependence of energy on stored charge for the 
classical metallic Coulomb Blockade device (blue dotted 
line) is changed to a discontinuous, and sometimes multi- 
valued, approximate parabola in the highly doped silicon 
single electron device (black solid line). Clearly the latter 
case will result in aperiodic Coulomb oscillations. More- 
over, the multiple values corresponding to particular elec- 
tron numbers result in the position of a Coulomb peak de- 
pending on which electron configuration is adopted (the 
energy barrier between configurations may prevent simul- 
taneous occupation). In this case, the Coulomb peak pat- 
tern can depend on the gate voltage history prior to the 
observation. 



IV. TRANSPORT UNDER MICROWAVE 
IRRADIATION 

In a typical experiment, the SET is biased with a fixed 
source-drain voltage and held at a fixed gate voltage, 
while the source-drain current is measured as the mi- 
crowave frequency is swept. This must be done with 
the SET set up as an impedance matched load (or with 
a matching impedance very close to the SET) or with 
the SET in a low quality factor (Q) cavity, which it- 
self is matched to the waveguide impedance. In the lat- 
ter case, the coupling efficiency is very greatly reduced 
and this approach has been used in this work. A typ- 
ical result is shown in Fig. |4j where the following fea- 
tures may be identified; large amplitude, low Q fluctua- 
tions (main plot), medium amplitude, mid Q fluctuations 
(right inset), and low amplitude, high Q resonances (left 
inset). The high Q features have been extensively stud- 
ied by Cresswell and coworkers, who have demonstrated 



T I I r 




Q I I 1 I 1 I I I 

0.80 1.00 1.20 1.40 1.60 



Frequency [GHz] 

FIG. 4: SET response to AC excitation. The large amplitude 
peaks have typically low Q value (few tens). Bias voltages: 
Vds=-|-0.7 mV, Vg=-2.77 V; radiation power: P=+15 dBm. 
Left inset: a high Q-value peak on an expanded frequency 
scale, with amplitude indicated by the double headed arrow. 
Q~10^. Right inset: resonances of different line-shape and 
medium Q closely spaced in frequency. Different data-set with 
respect to the main plot. Bias voltages: Vus^-l mV, Vg=- 
1 V; P= dBm. 



that the resonant frequency is affected by both the gate 
and the source-drain voltages, although the effect can be 
very small. Resonant features have previously been ob- 
served in two-dimensional electron gas systems (silicon 
MOSFETs and AlGaAs heterostructures) , where this be- 
haviour has been attributed to single particle excitation 
within the Fermi sea.l^In this previous work, a magnetic 
fleld is used to determine the energy difference between 
two levels, giving rise to a single resonance, resulting in 
only small numbers of resonances due to the requirements 
for a separately identiflable spin active system for each 
resonance. By contrast, the degenerately doped silicon 
SET exhibits large numbers of resonances without the 
need to apply a magnetic field. The large number of res- 
onances (many more than the expected number of free 
electrons in the island) suggests that their origin involves 
single particle excitation, with each able to give rise to a 
number of resonances. Furthermore, the photon energy 
corresponds to a small fraction of the thermal energy at 
the measurement temperature of 4K, so that the popu- 
lations in the two levels would be expected to be very 
nearly equal (of the order 0.49 and 0.51 for the upper 
and lower energy levels). 

The AC conduction behaviour of the hopping model, 
described earlier to explain the DC characteristics of de- 
generately doped silicon SETs, has been considered pre- 
viously.^"^ ^"^ All of these models consider statistically av- 
eraged behaviour over a large number of hops in order to 
predict a value for the conductance. However, in these 
SETs individual hops may contribute a significant change 
in the conductivity due to the small size of the device. 
The basic underlying mechanism is the change in hop- 



7 



ping rate between two localized sites, when exposed to 
electromagnetic radiation that is resonant with the en- 
ergy level difference between the sites. This changes the 
occupation levels at the two sites away from the values 
corresponding to thermal equilibrium. In an SET with 
localized sites, a change in just two sites may well give 
rise to a measurable change in the conductance, due to 
the very high sensitivity of a tunnel barrier transmission 
to wavefunction overlap. Any change in occupancy of 
localized sites close to the edge of a tunnel barrier will 
have an exponential effect on the transmission. If these 
changes in occupancy occur at sites further away, then 
the tunnel barrier is affected more indirectly through the 
polarisation of intermediate sites, which reduces the over- 
all change in transmission. 

Consider two sites separated by a distance that is sig- 
nificantly larger than the wavefunction overlap, and by 
a tunnel barrier that reduces the thermally driven tun- 
nelling rate between the two sites to a very low level, see 
left inset of Fig. [5j Under these circumstances, electro- 
magnetic radiation that is resonant with the separation 
between these energy levels will cause the electron in the 
occupied level to tunnel to and from the unoccupied level, 
at a rate (the Rabi rate) determined by the intensity of 
the radiation. If the intensity is very high, then the elec- 
tron will perform these spatial Rabi oscillations, as dis- 
cussed by Stafford and Wingreen.^^ Such oscillations can 
proceed coherently until there is an interaction with the 
environment that results in an exchange of energy and/or 
causes the electron to tunnel to a different site; such in- 
teractions determine the lifetime of these oscillations. 

The changes in site occupancy resulting from these spa- 
tial Rabi oscillations causes a change in the SET current 
by modifying the transmission properties of the tunnel 
barriers. This modification is a result of the change in 
electrostatic potential resulting from the electron dis- 
placement in the oscillation. The magnitude of this 
change depends on the separation between the sites, the 
distances and directions to the tunnel barriers and the 
Rabi rate. The electron displacement due to the Rabi 
oscillation gives rise to a time dependent polarisation, 
which gives rise to the change in potential at the tunnel 
barrier. However, this electric field is screened by the po- 
larisability of any intermediate structures, such as other 
tunnelling charges. If the spatial Rabi oscillation is slow, 
then screening can be very efficient and the affect on the 
tunnel barrier transmission is correspondingly small. As 
the Rabi rate is increased, the screening is progressively 
overcome and the tunnel barrier transmission is modi- 
fied. 

A single electron undergoing a spatial Rabi oscilla- 
tion will cause a monodirectional resonance (left inset 
of Fig. |4]). However, there is also a natural tendency for 
such systems to couple. The amplitude response of a sys- 
tem of coupled oscillators is well knowrP^ and results in 
a resonance whenever the driving field corresponds to a 
mode frequency. The individual motions of the electrons 
participating in the coupled mode is more complicated 



than the sum behavior. In particular, the relative ampli- 
tude for each coupled electron differs on either side of a 
resonant frequency. In cases where electrons acting inde- 
pendently would give rise to monodirectional resonances 
in opposite directions, the coupled motion would give rise 
to a differential behavior, as observed in Fig. |4] (right in- 
set). 

The width of each resonance directly refiects the life- 
time of the oscillation, which can vary over a very wide 
range. Rabi oscillations are usually limited to a max- 
imum rate that is very much slower than the resonant 
frequency of the excitation (say by a factor of 100). The 
lowest frequency at which high quality factor resonances 
have been observed in these experiments is of the order 
IGHz, so that the Rabi rate probably does not exceed a 
few tens of MHz. In addition, if the Rabi rate is slower 
than the rate at which the electron interacts with the 
environment, then complete oscillations cannot be made 
and there is no resulting polarisation to change the SET 
current. 

As briefiy mentioned above, the energy difference be- 
tween levels undergoing spatial Rabi oscillations, given 
by the resonant frequency, is much smaller than ksT 
at 4K, so that the Boltzmann function predicts a very 
small difference between the occupancies for these lev- 
els at thermal equilibrium. The spatial Rabi oscilla- 
tions, due to continuous wave radiation, can only change 
the occupancy levels towards the limiting case of 50% 
(right inset of Fig.[5|, so if these levels were already close 
to this figure there would be little change in polarisa- 
tion between the localized sites and so little change in 
device current. However, it is important to point out 
that for low dimensional structures at cryogenic tem- 
perature the electron-phonon interactions are strongly 
suppressed. I^ZISIl In particular, Tilke and coworkerS^Sihave 
found that for P-doped silicon nano-wires the electron- 
phonon relaxation time can be as high as few fis. Single- 
shot measurement^^ni performed on the same device sys- 
tem investigated here have shown that electron relax- 
ation time due to microwave excitation has a comparable 
timescale. Therefore, we can assume that the thermal- 
isation processes in the dot occur on a longer timescale 
than the photon-induced charge transfer does. This ac- 
counts for the observation of radiation effects at energy 
much lower than k^T. 

Finally, we turn to explain the effect of a variable ra- 
diation power on the resonant current. The Rabi rate, 
flji, increases as the square root of the microwave power 
under resonant excitation conditions. The magnitude of 
the device current difference is thus expected to change 
correspondingly. As pointed out earlier, the Rabi rate 
must exceed the rate of interaction with environment for 
any net charge displacement, so that in practice there is a 
threshold in microwave power that must be exceeded for 
any resonance to appear. The threshold condition occurs 
when Ti w 1 /fi_R where Ti is the relaxation time for the 
resonance. The occupancy of the normally unoccupied 
site increases with microwave power, saturating at 0.5 
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for very high powers. For intermediate power levels, the 
occupancy oscillates due to beating effects between the 
Rabi rate and Ti (see right inset of Fig. [s]) . The average 
occupation probability of the excited state is given by 



1 



Ti 



P= 7^ \ Pe{t)dt 

^1 Jo 



(9) 



where Pe(t) = sin'^{flfit/2) is the instantaneous occu- 
pancy. In order to evaluate the current dependence on ra- 
diation power (P) , we assume that each of the two states 
corresponds to currents Ig and /g. The instantaneous 
current is given by the weighted sum of the currents, 



I{t)=Igil-pe{t))+hpeit) 

and the average current is 



/ = mdt 

^1 Jo 



(10) 



(11) 



Since we are mainly interested in the current variation, 
we can set lg=0 which would result in 

AI = leP = le 7^ 

(12) 

In order to model the rapid suppression of any excita- 
tion in case Ti <C we introduce an exponential 
prefactor as follows 



A/ = 



1 



I _j_ p-a{Tifln-k) 

TiVLr/2 - cos{TiVlR/2)sin{TinRl2) 



(13) 



Tl^R 

with a, k fitting constants. As stated earlier, the Rabi 
rate is proportional to the square root of the radiation 
power and can be expressed as 



VLr = y/cP + 52 



(14) 



being c the constant of proportionality and 5 the detun- 
ing factor. Fig. [5] shows a comparison of the modelled 
excess current as a function of the radiation power with 
experimental data taken from a resonance at 1.9200 GHz 
and Q ~ 2000. As the power is swept, a non-monotonic 
behavior for the resonant current is observed for the vast 
majority of the radiation-induced peaks and the model 
successfully replicates its main features. However, at the 
lowest (highest) power levels the model response appears 
to be slightly slower (faster) than the real data response. 
This may be because the model does not include the ef- 
fects of charge screening and rectification which may play 
a role at low and high power regimes, respectively. 



V. CONCLUSION 

The very high doping density does not ensure metallic 
behavior in the tunnel barrier and the island of a sili- 
con SET. Increasing the doping density cannot overcome 



the sidewalls depletion causing the non-metallic behav- 
ior due to a consequential increase in sidewall trap sites. 
Transport through the tunnel barriers and island is domi- 
nated by hopping through a limited number of sites. Re- 
arrangement of a fixed number of electrons between a 
larger number of sites is seen to lead to aperiodic gate 
oscillations and varying peak heights. Individual trans- 
fers can be promoted by resonant microwave excitation, 
resulting in a spatial Rabi oscillation. Such transfers 
are detected through their effect on the DC transport. 
The equilibrium population and lifetime of these states 
is strongly infiuenced by the energetics of localized elec- 
trons. These effects result in each excitation correspond- 
ing to a different energy, in contrast to Zeeman splitting. 
These differences could be used to separately address ex- 
citations/electrons without the need for additional gates, 
greatly simplifying the problem of realising a multi-qubit 
quantum computer. Indeed, each radiation-induced res- 
onance could potentially be considered as a qubit, char- 
acterised by the resonant feature shape, frequency and 
width. The quantum system initial state would be de- 
fined by the set(s) of occupied donor sites that minimises 
the free energy in the system. Manipulation and gate 
operation could be based on well established NMR pro- 
tocols. 



< 




FIG. 5: Dotted line: experimental current dependence 
at resonance (Frequency =1.9200 GHz) for varying radiation 
power. A/_Ds is the excess current with respect to DC con- 
dition. Bias conditions: Vg=-2.0 V, Vi5s=-1.30 mV. Solid 
line: calculated current using eq. (n3| with /e=9.65xlO-''A, 

c=7.33xl0^^ ^^^1^^> 5=4.96x10^ rad/s, fc=4.98, a=104.04, 
Ti = lQ~^s. Left inset: schematic energy level diagram show- 
ing two donor locations and the occupied (filled circle) and 
unoccupied (empty circle) energy levels. Right inset: average 
probability of occupation of an excited state under resonant 
excitation against the ratio between the relaxation time and 
the Rabi period (eq.||9|). 



9 



Acknowledgments with D.E. Kmclnitskii. The authors would hke to thank 

Stephen Gibhn for helpful discussions and careful reading 
AR acknowledges financial support from EPSRC and of ^^e manuscript. 
NPL. DGH would like to acknowledge useful discussions 



^ N. Zimmerman and M. Keller, Meas. Sci. Technol. 14, 1237 
(2003). 

^ A. Steane, Rep. Prog. Phys. 61, 117 (1998). 
^ T. A. Fulton and G. J. Dolan, Phys. Rev. Lett. 59, 109 
(1987). 

"* L. Ji, P. Dresselhaus, S. Han, K. Lin, W. Zheng, and 

J. Lukens, J. Vac. Sci. Technol. B 12, 3619 (1994). 
^ R. Smith and H. Ahmed, J. Appl. Phys. 81, 2699 (1997). 
® A. Fujiwara and Y. Takahaslii, Nature 410, 560 (2001). 

N. Stone and H. Ahmed, Appl. Phys. Lett. 73, 2134 (1998). 
® J. Gorman, D. G. Hasko, and D. A. Williams, Phys. Rev. 

Lett. 95, 090502 (2005). 
^ R. Augko, W. Ebcrhardt, C. Single, F. Prins, D. Wharam, 

and D. Kern, Appl. Phys. Lett. 76, 2065 (2000). 
^° A. Tilke, R. H. Blick, H. Lorenz, J. P. Kotthaus, and D. A. 

Wharam, Appl. Phys. Lett. 75, 3704 (1999). 
" C. Single, F. E. Prins, and D. P. Kern, Appl. Phys. Lett. 

78, 1421 (2001). 

T. Altebaeumer and H. Ahmed, Jap. J. Appl. Phys. 42, 
414 (2003). 

C. Single, R. Augkc, F. E. Prins, D. Wharam, and D. P. 
Kern, Semicond. Sci. Technol. 14, 1165 (1999). 
^'^ M. Manscher, M. Savolainen, and J. Mygind, Trans. Appl. 
Supercond. 13, 1107 (2003). 

L. CreswcU, D. Hasko, and D. Williams, J. Appl. Phys. 
105, 104508 (2009). 
^® L. Creswell, D. Hasko, and D. Williams, Micro. Eng. 84, 
1614 (2007). 



^"^ T. Fcrrus, D. Hasko, Q. Morrissoy, S. Burge, E. Freeman, 

M. French, A. Lam, L. Creswell, R. Colhcr, D. Williams, 

ct al., J. Appl. Phys. 106, 33705 (2009). 

S. Vitkavage and E. Irene, J. Appl. Phys. 68, 5262 (1990). 
^® B. Shklovskii and A. Efros, Electronic properties of doped 

semiconductors ( Springer- Verlag, Berlin, 1984). 
^° D. V. Averin and Y. V. Nazarov, Phys. Rev. Lett. 65, 2446 

(1990). 

2^ G. P. Tribcris, J. of Phys. C: Solid State Physics 21, L821 
(1988). 

T. Fujisawa, Y. Tokura, and Y. Hirayama, Phys. Rev. B 
63, 081304 (2001). 
23 M. Migliuolo and T. Castner, Phys. Rev. B 38, 11593 

(1988). 

E. Ritz and M. Dressel, J. Appl. Phys. 103, 084902 (2008). 
'^^ C. Stafford and N. Wingreen, Phys. Rev. Lett. 76, 1916 
(1996). 

2® A. French, Vibrations and waves (Stanley Thornes, 1971). 

L. Pescini, A. Tilke, R. Blick, H. Lorenz, J. Kotthaus, 

W. Ebcrhardt, and D. Kern, Nanotcch. 10, 418 (1999). 

O. Bourgeois, T. Fournicr, and J. Chaussy, J. Appl. Phys. 

101, 016104 (2007). 
2^ A. Tilke, L. Pescini, A. Erbe, H. Lorenz, and R. BUck, 

Nanotech. 13, 491 (2002). 
3° D. Hasko, T. Fcrrus, Q. Morrisscy, S. Burge, E. Freeman, 

M. French, A. Lam, L. Creswell, R. Collier, D. Williams, 

et al., Appl. Phys. Lett. 93, 192116 (2008). 



